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ABSTRACT 



The NAVSTAR Global Positioning System is a worldwide, 
all-weather, satellite positioning system capable of high 
accuracy real-time position determinations. The Applied 
Research Laboratories ( ARL: UT) , in conjunction with the 
Naval Surface Weapons Center (NSWC/DL) and other government 
agencies, conducted geodetic field tests of a government 
sponsored prototype receiver, the 114 1 0 0/GEOSTAR , in March 
1984. Data were acquired from four satellites by three 
receivers, with antennas located at known stations, over 
approximately six-hour periods each day. 

Based on the real-time solutions acquired on March 1, 8, 
and 9, 1984, absolute point position determinations have an 
average discrepancy of 7.5 meters with a one sigma repeat- 
ability of less than a meter. Relative positions were 
determined to an average relative accuracy of 1:9,700 for 
distances of 14 to 26 kilometers. 
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I. INTRODUCTION 



The Naval Surface Weapons Center (NSKC) , with sponsor- 
ship by the Defense Mapping Agency (DMA) , has developed a 
system for geodetic positioning with the use of the NAVSTAR 
Global Positioning System (GPS) [Ref. 1]. A government 
sponsored prototype receiver, the II 4 100/G20STAR , was devel- 
oped by Texas Instruments through contract with Applied 
Research Laboratories, The University of Texas at Austin 
(ARL : DT) [Bef. 2]. The receiver was tested by ARL:UT in 
conjunction with NSWC and other government agencies over 
short and medium base lines in March 1984. The Defense 
Mapping Agency Hydrographic/Topographic Center (DMAHTC) nas 
expressed interest in independent analyses of the data from 
these tests [Ref. 3]. 

Absolute point positioning, as it is used in this 
thesis, is the determination of a position with the use of 
one GPS receiver. Relative positioning is the determination 
of two positions relative to each other with the use of two 
receivers. The basic premise of real-time relative posi- 
tioning is that GPS observation errors are to a large extent 
systematic in the same general area at the same time. If 
such is the case, relative positioning by two receivers can 
locate positions relative to each other more accurately than 
absolute positioning alone. 

Accuracy standards for horizontal control require 
certain relative accuracies between directly connected adja- 
cent points. Third-order surveys are used for local hori- 
zontal control in hydrographic surveying and other local 
projects. Third-order, Class II standards require relative 
accuracies of at least 1:5,000 between directly connected 
adjacent points (e.g., 5 meters relative accuracy between 
stations separated by 25 kilometers). [Bef. 4] 
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Geodetic positioning with GPS has advantages over 
conventional surveying methods. Visibility between sites 
and fair weather conditions are not required. Real-time 
positioning has a special application in the case where the 
position of a station is needed in a short amount of time 
and before the data can be post-processed. An analysis of 
the real-time solutions of the GEOSTAR short and medium base 
line tests will give an estimate as to what degree of accu- 
racy the system can achieve in real-time. 
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II. BACKGROUND 



A. THE GLOBAL POSITIONING SYSTEM 

The Global Positioning System is a worldwide, all- 
weather, satellite positioning system. GPS can be thought 
of as consisting of three components; the Space Segment, the 
Control Segment, and the User Segment. The Space Segment 
will consist of eighteen satellites (with three spares) in 
six orbit planes when fully operational in late 1988 or 
early 1989 [Ref. 5]. The satellites will orbit the Earth in 
twelve-hour periods at altitudes of approximately 20,000 
kilometers. This constellation will provide four satellites 
in view at one time at elevations of greater than 20° above 
the horizon. [Ref. 6] The satellites' ground paths will be 
fixed. The time at which a particular satellite will be in 
view at a certain position will be about four minutes 
earlier each day due to the differences between solar and 
sidereal time. [Ref. 7] 

The GPS satellites have on orbit weight of 950 pounds. 
Power is supplied by solar panels and nickel-cadmium 
batteries. The satellites are equipped with atomic clocks 
with stabilities on the order of 1 part in 10* 3 . [Ref. 8] 
Perturbations in the satellite orbits are caused by the 
forces of attraction of astronomical bodies, gravitational 
irregularities of the earth, and solar radiation pressure. 
These perturbations are compensated for by the Control 
Segment and freguency standard offsets. [Ref. 7] 

The Control Segment consists of four monitor stations 
located at Alaska, Hawaii, Guam, and Vandenburg AFB; and a 
Master Control Station (HCS) and an Upload Station at 
Vandenburg AFB [Ref. 9]. The monitor stations track the 



satellites then send the data to the MCS for processing to 
determine updated predictions of satellite and clock 
behavior. The updated predictions are transmitted to the 
satellites by the Upload Station. [Ref. 10] 

B. TRANSMISSION SIGNALS 

The GPS NAVSTAR satellites simultaneously transmit navi- 
gation information on 2 RF frequencies; the LI frequency at 
1575.4 MHz . and the L2 at 1227.6 MHz. Both frequencies are 
modulated with pseudo-random noise (PRN) codes. The Precise 
Positioning Service (PPS) , formerly called the "P-code", is 
modulated on both frequencies and the Standard Positioning 
Service (SPS) , formerly called the "C/A-code", is available 
only on the LI frequency. PRN codes, different for each 
satellite, serve the purpose of uniquely identifying the 
particular satellite ana provide a means of determining the 
signal transit time by measuring the phase shift necessary 
for the receiver to match codes. [Ref. 11 and 12] 

The SPS is a relatively short code of 1023 bits. The 
code repeats itself every millisecond at a bit rate of 1.023 
Mbps. Due to the short duration time of the SPS, the code 
is relatively easy to match and lock on to. The PPS has a 
7-day code duration and 10 times the bit rate, making the 
PPS code relatively difficult to match with the receiver’s 
internally generated code. The PPS offers the advantage of 
a finer resolution of the phase shift necessary to match 
codes, thus a greater measurement accuracy of signal transit 
time. [ Ref. 11] 

C. THE NAVIGATION MESSAGE 

The GPS Navigation Message is modulated on both the LI 
and L2 frequencies and contains information necessary to 
compute the user’s position. The information is in the form 
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of a 50 tps. data stream common to both tne PPS and SPS . 

The data frame is 1500 bits long and consists of 5 subframes 
of 6 seconds (300 bits) each. The Navigation Message 
provides information such as satellite e^hemerides, system 
time, satellite clock behavior, and transmitter status. 

[Ref. 13] 

Subframes are divided into ten 30-bit words. The first 
two words of each subframe contain system time and time 
synchronization information for transfer from the SPS to the 
PPS in the form of hand-over words (HOW) and telemetry words 
(TLM) . The remaining 8 words contain user navigation data 
generated by the Control Segment. [Ref. 13] 

Data 3lock 1 occupies words 2 through 10 in subframe 
one. It contains frequency standard corrections, an associ- 
ated age of data word, and ionospheric delay corrections for 
single frequency receivers. The clock correction parameters 
are used to determine the satellite’s clock offset from GPS 
time by the following equations: 



t 



t 

sv 



At 



sv 



( 2 . 1 ) 



At = a + a, (t - t ) ♦ a„ (t - t ) 2 (2.2) 

where t is the GPS time in seconds measured from the begin- 
ning of the GPS week (approximately Saturday/Sunday midnight 
GMT), t is the satellite’s PRN code phase time at the tine 
of message transmission, t QC is Data Slock 1 reference time 
(seconds), and a Q , a 1 , and a 2 are Data Block 1 clock correc- 
tion parameters (coefficients of a second order polynomial) . 
Note that t may be approximated by t in equation 2.2. 

[Ref. 13] 

The age of data word (AODC) is the time between Data 
Block 1 reference time and the time since the last measure- 
ment used in computing the clock correction parameters 
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[Ref. 13]. Test data in this thesis were obtained using 
both the LI and L2 frequencies, therefore, the ionospheric 
delay parameters for single frequency receivers contained in 
Data Block 1 are not used (see Van Dierendonck, et al, 

1978) . 

Data Block 2 contains ephemeris information in subframes 
2 and 3. This data includes associated age of data words 
(AODE) which represent the time difference between Data 
Elock 2 reference time and the time of the last measurement 
used in the estimation of the ephemeris parameters. The 
ephemeris data includes conventional Keplerian-type parame- 
ters, secular drift terms, and harmonic coefficients used to 
correct for perturbations in the Kaplerian orbit. [Ref. 13] 
The satellite's position is calculated by Van Dierendonck ' s 
algorithm [Ref. 13] with the addition of a secular drift 
term to correct for perturbations in the inclination of the 
satellite's orbit [Ref. 14]. 

The Message Block appears in the third through the tenth 
30-bit words of subframe 4 of the Navigation Message. This 
space is reserved for alphanumeric messages provided to the 
user by the Control Segment. Space exists for 23 eight-bit 
ASCII characters. [Ref. 13] 

Almanac information appears in Data Block 3 in the fifth 
subframe. Data Block 3 consists of 25 subframes allowing 
for almanacs of up to 24 satellites (with one dummy 
subframe) . Twenty-five navigation messages must be received 
to acquire the complete data set of Data Block 3. Extra 
subframes can be used to repeat the almanacs for certain 
satellites. Satellite almanacs provide the user with less 
precise position and clock correction information. Almanacs 
aid the user in satellite selection and simplify the acqui- 
sition of signals. Almanacs include ephemeris parameters, 
clock correction parameters, satellite identification, and 
satellite health terms. [Ref. 13] 
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RECEIVER OPERATIOH 



The GEOSTAR GPS receiver is a government sponsored 
version of the TI4100 commercial receiver (Fig. 2-1). The 
GEOSTAR is a single channel, multiplex receiver capable of 
tracking four satellites on the LI and L2 PPS. The control 
and display unit (CDU) is a hand-held input/output device 
with two 16-character alphanumeric displays. The antenna 
assembly consists of a conical omnidirectional antenna with 
a preamp that allows remote operation of the receiver. Two 
100-foot cables connect the antenna assembly to the 
receiver. One cable is used to input DC power to the preamp 
and output the amplified EF signals. The other cable is 
used for a built-in test signal injection from the receiver. 
[Ref. 15] 

A typical GEOSTAR geodetic surveying unit consists of a 
2-man team, a station wagon or van, and the GSOSTAR field 
set. The field set contains the receiver, CDU, the antenna 
assembly, a cassette tape recorder, a tripod with a leveling 
bracket, two 12-volt automobile batteries, an EFRATOH FRK-H 
rubidium oscillator, a weather station with sensors, cables, 
transportation cases, tools, manuals, and spare parts. The 
field set is easily contained and transported in the rear of 
the vehicle. 

The antenna assembly is mounted on a standard surveyor’s 
tripod with a leveling bracket then leveled and plumbed over 
the mark to be surveyed. A reference mark on the base of 
the antenna assembly is oriented to north using a hand 
compass. The height of the L2 phase center mark at the base 
of the antenna assembly above the survey point is measured 
and recorded. The antenna assembly, the power supply, the 
cassette recorder, the weather station, and the rubidium 
oscillator are then connected to the receiver with the 
appropriate cables. Scratch tapes are inserted in the 
cassette recorder for the power-up procedure. 
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Figure 2.1 The TI4 100/3EOSTAR Receiver 
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The system is turned on and the receiver and cassette 
recorder self-tests are initiated. ilhen the self-tests are 
successfully completed, the GESAE bootstrap tape is inserted 
in Drive 1. The program is loaded and blank tapes are 
inserted in the two tape drives. The operator is then 
prompted by the CDU for input information including an esti- 
mate of the user's position, weather data, satellites to be 
tracked, ana Doppler aiding values for each satellite. The 
system then operates unattended except for the replacement 
of cassette tapes as required and the periodic manual input 
of weather data (if the automatic weather station is not 
used) . 



E. POSITION COHPOTATION 
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F. GEODETIC DATUMS 



The GEOSTAR real-time solutions are referenced in DoD World 
Geodetic System 1972 (WGS-72) Cartesian coordinates. WGS-72 

is an earth-centered, earth-fixed, glonal, geodetic datum 
adopted by DMA in March 1974 [Ref. 18]. The reference 
surface is an ellipsoid that test approximates the shape of 
the earth’s geoid. The ellipsoid is normally defined by the 
length of the semimajor axis and the flattening. 

Latitude ranges between -90° and 90° with north posi- 
tive. Longitude is measured positive east from 0° to just 
less than 360° from the WGS-72 reference meridian at 0".54 
east of the Greenwich Meridian. The Cartesian X-axis is 
formed by the intersection of the reference meridian plane 
and the plane of the equator. The Y-axis is 90° east of the 
X-axis in the plane of the equator. The Z-axis is in the 
direction of the earth's mean rotation axis based on the 
average terrestrial pole of 1900-05. [Ref. 19] 

The North American Datum 1927 (NAD-27) is a regional 
datum based on the Clarke 1866 Ellipsoid with the origin at 
Meades' Ranch, Kansas. The Clarke 1366 Ellipsoid minimizes 
tne deflection of the vertical and geoidal undulations in 
the continental United States. The National Ocean Survey 
(NQS) , the National Geodetic Survey (NGS) , and the United 
States Geological Survey (USGS) use NAD-27 coordinates in 
most survey operations. The transformation from NAD-27 to 
WGS-72 reference systems creates offsets of the order of 5 
meters [Ref. 2]. It is assumed that this error is removed 
when the reverse transformation is made; but the error 
source will still exist in the analysis of the WGS-72 point 
position discrepancies. 
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III. ANALYSIS METHODS 



A. EXPERIMENT DESIGN 

Relative position field testing of the GEOSIAR Receiver 
was performed in central Texas in February and March 1984. 
Three receivers, each with an antenna located over a 
different known position, were operated continuously over 
approximately the same period of satellite visibility. Four 
satellites, with PEN codes 6, 8, 9, and 11, were tracked in 

these tests. The period of data gathering was coordinated 
with the time period in which all four satellites were in 
view at the same time. The real-time positions of the 
antennas were recorded every 12 seconds over approximately a 
six-hour period each day. The data were recorded on 
cassette tape then transferred to 9-track tape by AEL:UT. 

B. GPS USER TAPES 

The recorded data and the associated position coordi- 
nates used in the analysis were obtained from NSWC/DL. The 
GPS User Tapes are a standardized format, 9-track magnetic 
tape. The tapes are unlabelled, ASCII coded, 1600 bpi, with 
80 characters per physical record and 1 record per block. 

The first file on the tapes consists of header comment 
records. This file assists the user in identifying the tape 
and includes information such as the date and site of the 
data acquisition. 

Data files begin with an 80-character control record 
(Data Item List) which identifies the types of data that 
follow. The first location of all records is reserved for a 
control character. The Data Item List control record has an 
asterisk 1 in the first position followed by an integer in 
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the next three positions which specifies the nuicer of times 
to repeat the sequence of the Data Item List. The remaining 
76 characters contain up to 19, 4-digit data type identi- 
fiers which specify the type of data that is to follow the 
control record. Control records are read with a 
(A1, 13, 1914) format. 

Data information are one or more records in length and 
divided into groupings called data items (Tab. I) 

£Eef. 19]. 



TABLE I 

Data Types of GPS Jser Tapes 



Item Ide nti fier Data Type 



000 1-0999 
1000-1999 
2000-2999 
3000-3999 
4000-4999 
5000-5999 
6000-6999 
7000-7999 
8000-8999 



Exchange conversion program information 

Constant equipment related data 

Campaign identifying data 

Infrequent measurement data 

Space vehicle related data 

Weather data 

Real-time solution data 

Measurement data 

Post analysis data 



The data type identifiers of the Data Item List specify the 
order and the format of how the data is to be read. 

Data used in the analysis were read from the user 
exchange tapes using a slight modification of the sample 
input program contained in reference 19 (App. A) . 

Real-time solution information (data type 6010) was read 
from the tapes and written to mass storage and retained. 
Other information such as measurement data types (data type 
920), real-time solution configuration (data type 1150), and 
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anter,na location (data type 2210) were read manually from 
the first page of the printouts of the tapes. 

C. COORDINATE TRANSFORMATIONS 

The coordinates of the stations used in the field 
testing were given in NAD-27 geodetic coordinates. The 
real-time solution output of the GEOSTAE Receiver are in 
WGS-72 Cartesian coordinates. To make comparisons of the 
given positions to the observed positions, the given posi- 
tions were transformed to WGS-72 Cartesian coordinates. The 
transformation was done using the Abridged Molondensky 
formulas contained in the NGS's conversion program (D035-41) 
for the HP-41CV programmable calculator. The program 
prompts the user for the values of the origin shift between 
coordinate systems. The X, Y, and Z origin shifts for 
transforming from NAD-27 to WGS-72 are -22, 157, and 176 
meters, respectively [Ref. 20]. These values differ from 
the origin shift values supplied with the transformation 
program by NGS. The NGS values are -12.547, 155.804, and 

175.369, for the X, Y, and Z origin shifts, respectively. 

The program generates semimajor axis and inverse flattening 
parameters of both reference systems (Tab. II) . 



TABLE II 

Transformation Program Parameters 



El lips oid 



a (me ters ) 



1/f 



Clarke 1866 



6,378,206.4 

6,378,135.0 



294. 9786952 



WGS 72 



298.2600000 



D. SATELLITE UPDATES 

The ephemeras elements and clock model updates occurred 
at identical times according to the data recorded in 
field tests. An approximation of the update times is calcu- 
lated by subtracting the Age of Data (AODE) from the ephem- 
eris reference time (t ) , i.e., 

update time = t - AODE (3.1) 

oe 

The update times changed by -496/ 1052, or 1552 seconds when 
the computation of equation 3. 1 was performed for each 
satellite with subsequent occurrences of the broadcast 
ephemeris on the recording. The change is due to packing of 
bits in the transmissions of the AODE values. An approxima- 
tion of the update times is made, within a few minutes, by 
using the first AODE values that cnange to smaller numbers 
(on the order of 1 to 2 hours) from previously recorded 
large values (on the order of 24 hours) . [Hef. 14] 

E. AVEEAGING METHODS 

Relative position determinations were made using real- 
time solutions from two different stations at the same GPS 
time tags. Computer programs were written to read the real- 
time solutions from mass storage and compute the discrepan- 
cies between the observed and the given WGS-72 coordinates 
of pairs of stations at common GPS times. The beginning and 
ending times of data acquisition varied among the three 
receivers operated during the same time period; thus, the 
effect of comparing solutions at the same time tags resulted 
in the elimination of some of the real-time solutions from 
the analysis. In effect, the tape with the earlier starting 
time was spooled forward to correspond to the starting time 
of the first real-time solution recorded by the receiver 
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that began operation at a slightly later time. Similarly, 
the time of the last solution common to both receivers is 
the last solution recorded by the receiver that ended opera- 
tion first. 

A mean position was determined in two different ways to 
show the effect of the ephemeris and clock updates on the 
real-time solutions and to eliminate real-time solutions 
with large discrepancies from the averaging. In general, 
relatively large discrepancies are evident in roughly the 
first hour of receiver operation. The real-time solution 
discrepancies of pairs of stations were averaged beginning 
with the values one hour after and two hours after the first 
solution time common to both receivers. In these tests, 
beginning the averaging one hour after the initial solutions 
corresponds to averaging the solutions that were recorded a 
short time before or after the update times and through the 
remainder of the recording session. Beginning the averaging 
two hours after the initial solutions corresponds to using 
only those values that were recorded much after the updates 
and through the remainder of the session. Real-time solu- 
tions are treated as independent observations in the statis- 
tical analysis. 

F. COMPUTATIONS 

The distances between stations (d) , in the WGS-72 refer- 
ence system, are computed by equation 3.2 where X ^ , Y]_, and 
Z 1 are the WGS-72 Cartesian coordinates of one of the 
stations, and X 0 , Y 0 , and Z 0 are the coordinates of the 
other station of the side. The side-lengths are calculated 



d 




(3.2) 
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foe both the given coordinates and the mean values of the 
GEOSTAR real-time solution coordinates. The discrepancies 
between the two are computed by subtracting the given 
distances from the observed distances. Discrepancies are 
used in the averaging to save computer time and to reduce 
round-off error. 

The absolute value of the point position discrepancies 
are computed by equation 3.3 where the subscripts o and g 
represent the observed and given values, respectively. 

6 = T/(X o - X g >* ♦ d 0 - I,)* * <z„ - Z,)= (3.3) 

The standard error of a single observation was calcu- 
lated as a measure of the precision of the real-time solu- 
tions. The means and standard errors of a single, 12-second 
real-time solution were calculated for the coordinate 
discrepancies and the point position discrepancies. 
Comparisons are made for the two averaging methods (Sect. 

Ill E) . The mean WGS-72 Cartesian coordinate discrepancies 
were used to compute the final observed positions. The 
observed positions are computed by subtracting the discrep- 
ancies from the given coordinates. 

G. GEODETIC COHPARISONS 

Geodetic accuracy standards for horizontal control 
require certain relative accuracies between directly 
connected adjacent points. The distances between points are 
computed along the ellipsoid - the geodesic. To compare the 
accuracies of the observed positions with the standards for 
geodetic control, the WGS-72 Cartesian coordinates of the 
observed positions are transformed to NAD-27 geodetic coor- 
dinates. The transformations were done using a modification 
of NGS’s D035 program for the HP-41 C V programmable 
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calculator. The modified version uses the same parameters 
as those contained in the original program (Tab. II). 

The geodetic inverse computation is used to calculate 
the geodetic distance and azimuths between two known points. 
The inverses between given positions and between observed 
positions were computed and compared. 

To simulate the positioning of an unknown station, one 
station of a side is held fixed and the other station is 
located using the direct computation. The coordinates of 
the known position and the distance and azimuth to the 
unknown position are the input variables and the coordinates 
of the unknown station are computed. The distance and 
azimuth calculated by the inverse computation of the 
observed point positions are used in this calculation. The 
direct and inverse computations were done on the IBM 3033 
based on Puissant’s Coast and Geodetic Survey Formula 
[Ref. 21 ]. 

The standard errors of the discrepancies of the geodetic 
position determinations for entire recording sessions are 
calculated as an approximation of the repeatability. To do 
so, the two point position discrepancies of each station on 
days 61 and 63 must first be averaged to determine a single 
value for the station. 
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IV. RESOLTS AND DISCUSSION 



A. GIVEN POSITIONS AND DISTANCES 

Data in this report are from relative position tests of 
March 1, 8, and 9, 1984. Three receivers were operated 
simultaneously over a period of approximately six hours. 

The antennas of the March 1st tests were located over a 
second-order position and two of its reference marks. The 
triangle has side-lengths of approximately 30 meters. The 
tests of March 8th and 9th were conducted at first-order 
survey points separated by distances of 14 to 25 kilometers. 
The data of the tests at one of the positions on March 9th 
was not available at the time of this writing. The data 
gathered on March 1st will be defined as the short base line 
tests and the 14 to 25 kilometer testing will be termed 
medium base line tests. Tables III and IV give the names 
and the NAD-27 coordinates of the stations and Table V gives 
the distances and azimuths in each triangle. Distances and 
azimuths were calculated using the inverse computation 
(Sect. Ill G) . 

The NAD-27 coordinates of the given positions were 
transformed to WGS-72 Cartesian coordinates using a NGS 
transformation program (D035-41). The ellipsoidal heights 
are corrected for the heights of the antennas above the 
station marks prior to the transformations. The heights of 
the antennas above the stations (data type 2110) were read 
manually from the first page of the printouts of the tapes. 
The distances between the points in N3S-72 Cartesian coordi- 
nates were calculated using equation 3.2 (Tab. VI). 

The computed side distances for the FGS-72 Cartesian 
coordinate system are in all cases greater than or equal to 
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In the medium base line 



the computed geodetic distances, 
tests , the differences were found to increase by . 0035 



TABLE III 
Station Names 



J.D. 


Number 


Name 




Abbreviation 


61 


85019 


ARL 10563 


DMA 


JMR 




35020 


ARL 10563 


RM 1 DMA 


PM 1 




35021 


ARL 10563 


RM 2 DMA 


RM 2 


68,69 


35022 


BP. ZERO 




BZERG 




35027 


PILOT KNOE 


1935 


PILOT 




* 85028 


BALD KNOB 


1935 


BNOB 



* Station not used in March 9th analysis 



TABLE 17 

NAD-27 Station Coordinates 



Station 


North Latitude 


West 


L2Q.3.i.tude 


Ht . (m.) 


JMR 


300 


22 ’ 


56". 74285 


970 


43’ 


35". 97558 


235. 5 


RM 1 


30° 


22 ’ 


57". 54452 


970 


43’ 


36". 31041 


235. 7 


FM 2 


300 


22 ' 


57". 06234 


970 


43' 


35 ". 0 1009 


235.6 


BZERO 


300 


23’ 


16". 31238 


970 


43’ 


42'*. 05790 


238. 4 


PILOT 


300 


09’ 


30". 93387 


970 


42' 


28". 08922 


215. 9 


EN0B 


300 


20 ' 


07". 61000 


97° 


35' 


30". 84700 


204. 4 


percent 


of the 


len 


gth of the 


line. 


Differences 


of .911, 



.748, and .534 meters were computed for the side-lengths of 
approximately 25.5, 22.6, and 14.3 kilometers, respectively. 
These differences occur because the Cartesian distances are 
computed between each point; the geodetic distances are 
computed from the projections of the points on the ellipsoid 
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which does not account for the ellipsoidal heights- The 
distances between BZERO and PILOT are equal to the nearest 



TABLE V 

NAD-27 Geodetic Distances and Azimuths 



Side 




Distance (m. ) 


Azimuth 




JMR - 


RM1 


26.255 


3400 


05’ 


34". 72573 


J ME - 


RM2 


27. 59 1 


690 


06’ 


37" . 54898 


EMI - 


RE2 


37.759 


1130 


09' 


19". 14446 


BZERO 


- PILOT 


25, 492.277 


1750 


32’ 


48". 25379 


BZERO 


- BNOB 


14,346. 955 


1130 


51' 


26" . 80840 


BNOB - 


• PILOT 


22,556. 303 


209° 


40 ’ 


09". 53915 





TABLE VI 






WGS-72 Side-length 


Distances 


J. D. 


Side 


Distance (m. 


61 


JMR - R M 1 


26.255 




JMR - R M2 


27. 593 




RM 1 - RM2 


37.759 


69 


BZERO - PILOT ‘ 


25, 493. 183 




BZERO - BNOB 


14, 347.439 




BNOB - PILOT 


22, 557.051 


69 


BZERO - PILOT 


25, 493 .188 



millimeter for tests run on days 68 and 59 even 
heights of the antennas above the stations were 
both days. 



though the 
different on 
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B. SATELLITE UPDATE TIMES 



The ephemeris ana clock update tines occurred approxi- 
mately one hour after the start of data acquisition in these 
tests. There are approximately 25 to 30 recordings of the 
Lroadcast ephemeris (data type 4110) during the six-hour 
recording periods. The update times were computed (eqn. 

3.1) using the AODE values that changed from large numbers 
to much smaller numbers in subsequent recordings of the 
Lroadcast ephemeris (Tab. VII). Satellites with PSN codes 
9 and 11 were updated approximately one-half hour after the 
updates of codes 6 and 8 on days 61 and 69. According to 
the recorded data, all four satellites were updated at 
approximately the same time on day 68. 



TABLE VII 

Satellite Update Times (Seconds) 



J. D. 


PPN codes t ^ 

oe 


AODE 


- A0 


61 


6,8 


378,000 


6, 144 


371, 856 




9,11 


378,000 


4,096 


373, 904 


68 


6,3,9, 1 1 


378,000 


6,144 


371,856 


69 


C R 
•" / 


460,300 


4, 096 


456, 704 




9,11 


464,400 


6,144 


458,256 



C. BEHAVIOR OF THE REAL-TIME SOLUTIONS 

1 • WGS -72 Coordinate Discrepancies 

The real-time solutions were read from the GPS User 
Tapes and retained. The solutions of pairs of stations, 
with time tags of real-time solutions common to both, were 
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used in the analysis. The discrepancies between the 
observed coordinates and the given WGS-72 coordinates 
(observed - given) were computed and plotted over time 
(Figs. 4.1 - 4.8). The vertical arrows in the figures 
indicate the times of satellite updates (Tab. VII). The 
tick marks on the horizontal (time) axis represent one hour. 

In general, the X, Y, and Z coordinate discrepancies 
begin with large values in the initial solutions and then 
become smaller and more consistent after satellite updates 
and over time. There is a significant improvement in the 
discrepancies of the observed coordinates near the time of 
the update of PEN codes 6 and 8 at station JdR and near the 
update of codes 9 and 11 at station P.'ll on day 61 (Figs. 

4.1 and 4.2). The precision of the discrepancies in each 
coordinate increase at all three stations on day 61, a snort 
time after the update of codes 9 and 11. The discrepancies 
at stations BZEEO and PILOT on day 68 decrease approximately 
one hour after the initial solutions (Figs. 4.4 and 4.5) . 

2. WGS-72 Point Position D is crepancies 

The absolute values of the point position discrepan- 
cies of pairs of stations with common time tags were 
computed (eon. 3.3) for five-minute time segments and 
plotted with the standard errors over time (Figs. 4.9 - 
4.16). In general, the point position discrepancies and the 
standard errors decrease after satellite updates and over 
time. Except for the general trend, there is no obvious 
correlation of the accuracies (discrepancies) and the preci- 
sion of the five-minute time segments at all stations on the 
same day. 



3. WGS-72 Sid e-len gth Discrepancies 



relatively high accuracy side-length determinations 
are evident at certain times; however, the times of these 
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DISCREPANCY (METERS) 

- 70 . 0 - 60.0 - 50.0 - 40.0 - 30.0 - 20.0 - 10.0 0.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 



J.D. 61, 1984 



SHORT BASE LINE TESTS 




Figure 4.1 WGS-72 Coordinate Discrepancies 
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DISCREPANCY (METERS) 

- 70 . 0 - 60.0 - 50.0 - 40.0 - 30.0 - 20.0 - 10.0 0.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 



J.D. 61, 1984- 

SHORT BASE LINE TESTS 
SAMPLING INTERVAL = 5 MINUTES 




Figure 4.2 WGS-72 Coordinate Discrepancies 
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DISCREPANCY (METERS) 

- 70 . 0 - 60.0 - 50.0 - 40.0 - 30.0 - 20.0 - 10.0 0.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 



J.D. 61, 1984 

SHORT BASE LINE TESTS 
SAMPLING INTERVAL = 5 MINUTES 




Figure 4.3 WGS-72 Coordinate Discrepancies 
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DISCREPANCY (METERS) 

- 70 . 0 - 60.0 - 50.0 - 40.0 - 30.0 - 20.0 - 10.0 0.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 



J.D. 68, 1984 

MEDIUM BASE LINE TESTS 
SAMPLING INTERVAL = 5 MINUTES 




Figure 4.4 HGS-72 Coordinate Discrepancies 
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DISCREPANCY (METERS) 

- 70 . 0 - 60.0 - 50.0 - 40.0 - 30.0 - 20.0 - 10.0 0.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 



J.D. 68, 1984 



MEDIUM BASE LINE TESTS 




Figure 4.5 WGS-72 Coordinate Discrepancies 
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DISCREPANCY (METERS) 

- 70 . 0 - 60.0 - 50.0 - 40.0 - 30.0 - 20.0 - 10.0 0.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 



J.D. 68, 1984 

MEDIUM BASE LINE TESTS 
SAMPLING INTERVAL = 5 MINUTES 



i — 




Figure 4.6 WGS-72 Coordinate Discrepancies 
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DISCREPANCY (METERS) 

- 70 . 0 - 60.0 - 50.0 - 40.0 - 30.0 - 20.0 - 10.0 0.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 



J.D. 69, 1984 

MEDIUM BASE LINE TESTS 
SAMPLING INTERVAL = 5 MINUTES 







Figura 4.7 WGS-72 Coordinate Discrepancies 
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DISCREPANCY (METERS) 

- 70 . 0 - 60.0 - 50.0 - 40.0 - 30.0 - 20.0 - 10.0 0.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 



J.D. 69, 1984 

MEDIUM BASE LINE TESTS 
SAMPLING INTERVAL = 5 MINUTES 




Figure 4.8 HGS-72 Coordinate Discrepancies 
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I 



STATION JMR 
J.D. 61,1984 




Figure 4.9 Discrepancies and Standard Errors 
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STATION RM1 
J.D. 61,1984 




Figure 4.10 Discrepancies and Standard Errors 



1 



STATION RM2 
J.D. 61, 1984- 




Figure 4.11 Discrepancies and Standard Errors 
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STATION BZERO 
J.D. 68,1984 




Figure 4.12 Discrepancies and Standard Errors 
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STATION PILOT 
J.D. 68,1984 




Figure 4.13 Discrepancies and Standard Errors 



43 



STATION BNOB 
J.D. 68,1984 




Figure 4.14 Discrepancies and Standard Errors 
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STATION BZERO 
J.D. 69,1984 




Figure 4.15 Discrepancies and Standard Errors 
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METERS 



STATION PILOT 
J.D. 69,1984 




Figure 4.16 Discrepancies and Standard Errors 
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solutions are not consistent at all stations luring the same 
day (Figs. 4.17 - 4.23). The discrepancies in the lengths 
of sides JMR-RM2 and EM1-P.M2 are relatively small near time 
tag 382632; however, the side- length discrepancy of JMR-RM1 
at approximately the same time is relatively larger than 
other determinations of that side on day 61. The side- 
length discrepancies of JMR-RM2 and RM1-RM2 are -.044 and 
.037 meters at time tags 382632 and 382600, respectively. 

The discrepancy of JMR-RM1 at time tag 382600 is .955 
meters; whereas, the discrepancy is .476 meters at time tag 
386400. Therefore, the most accurate side-length determina- 
tions do not necessarily occur at approximately the same 
time for all three sides. 

The most accurate side-length determinations do not 
necessarily occur when the magnitudes of the point position 
discrepancy vectors of two stations are nearly the same 
because of differences in the directions of the discrepan- 
cies. The side-length discrepancies determined on day 68 
for the side BZEEO-BNOB illustrate this point (Fig. 4.21). 
The point position discrepancies for stations BZERO and 3N0B 
are 10.889 and 10.816 meters at time tag 377808 and the 
side-length discrepancy is -.151 meters. At time tag 
384408, the side-length discrepancy is -.036 meters and the 
magnitudes of the point position discrepancies are 9.295 and 
14.053 meters. The side-length discrepancy is smaller. at 
the later time tag even though the difference in the magni- 
tudes of the point position discrepancies are greater at 
that time. 

D. AVERAGING RESOLTS 

1 . ^termination of the Method of Averaging 

Two averaging schemes were used to arrive at mean 
position discr epancies (Sect. Ill E) . The means and 
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J.D. 61, 1984 



SHORT BASE LINE TESTS 




Figure 4. 17 Point Position and Side-length Discrepancies 
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J.D. 61, 1984 



SHORT BASE LINE TESTS 




Point Position and Side-length Discrepancies 
49 



Figure 4. 18 



ABSOLUTE VALUE OF DISCREPANCY (METERS) 



J.D. 61, 1984 

SHORT BASE LINE TESTS 



SAMPLING INTERVAL = 5 MINUTES 







Figure 4.19 Point Position and Side-length Discrepancies 
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J.D. 68, 1984- 

MEDIUM BASE LINE TESTS 




Figure 4.20 Point Position and Side-length Discrepancies 
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J.D. 68, 1984- 

MEDIUM BASE LINE TESTS 
SAMPLING INTERVAL = 5 MINUTES 




Figure 4.21 Point Position and Side-length Discrepancies 
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J.D. 68, 1984 

MEDIUM BASE LINE TESTS 




j 



Figure 4.22 Point Position and Side-length Discrepancies 
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J.D. 69, 1984 

MEDIUM BASE LINE TESTS 
SAMPLING INTERVAL = 5 MINUTES 







Figure 4.23 Point Position and Side-length Discrepancies 
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standard errors of discrepancies were computed for the X, Y, 
and Z coordinates and the absolute values of the total point 
position discrepancies (Tabs. VIII and IX). 

The standard errors of the point position discrepan- 
cies for the solutions that were averaged beginning two 
hours after the initial solutions are, in general, smaller 
than those that were determined using all solutions begin- 
ning one hour after the initial solutions. The average 
standard error of the absolute value of the point position 
discrepancies for all pairs of determinations is 2. 302 
meters using the former averaging method, compared to 1.422 
meters using the latter averaging method. Much of the 
difference between these two standard errors can be attrib- 
uted to the relatively large point position discrepancies 
computed for stations 3SEF.0 and PILOT on day 69 between the 
first and second hours after the initial solutions (Figs. 
4.15 and 4.16). Neglecting stations BZERO and PILOT on day 
69, the average of all standard errors of the point position 
discrepancies are 1.657 meters and 1.333 meters for the two 
different averaging methods. 

The elimination of the real-time solutions recorded 
between one hour and two hours after the start of data 
acquisition results in smaller standard errors of the point 
position discrepancies in the averaging process. The 
remainder of the analysis uses the mean coordinates deter- 
mined by the averaging method in which only the real-time 
solutions recorded two hours after the initial solutions are 
used (Tab. IX) . 

2 • WGS-72 R esul ts 

The Y coordinate discrepancies are larger than the 
discrepancies of X and Z coordinates at all stations with 
the exception of those at station BZERO on day 69. The oean 
observed X coordinates are consistently different from the 
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TABLE 

Means and Sigmas Beginning One 



Station 


X mean 
Sigma X 


Y mean 
Sigma 


JMR 


6. 677 
1.242 


-12. 696 
1. 558 


R M 1 


5. 178 
0.688 


-13. 512 
1. 927 


JMR 


7. 502 
3. 356 


-11.648 
4. 593 


RM 2 


6.205 
2. 574 


-12. 186 

5.639 


RM 1 


5. 162 
0.700 


-13. 623 
1. 900 


RM 2 


5. 588 

1. 165 


-13.548 
1. 375 


BZERO 


5. 806 
0. 725 


-7. 578 
3.818 


PILOT 


5. 653 
0.57 1 


-7. 877 
1. 107 


BZERO 


5. 972 
0. 987 


-7. 337 
1. 4 52 


B NOB 


6. 060 
1.625 


-9. 125 
5. 337 


PILOT 


5. 708 
0.59 1 


-7. 916 
1.178 


B NOB 


5.822 
1. 24 1 


-9. 777 
4. 527 


BZERO 


6.629 

1.379 


-1. 833 
4. 421 


PILOT 


6.613 
3. 050 


-1. 030 
11.727 



Dnits are meters. 



VIII 



Hour After 


Initial 


Solutions 


Z mean 
Sigma Z 


Point 

Position 


Sample 

Size 


3.196 

1.292 


14.834 
1 . 244 




5.371 

0.635 


15.495 

1.667 


151 8 


2.426 
3. 296 


15.372 

2.178 




4 . 789 
4.619 


16.274 

2.197 


1564 


5.405 

0.682 


15.600 

1.632 




5.794 

2.707 


16.019 

1.512 


1450 


-2. 623 
1.226 


10.006 

0.761 




1.090 

1.240 


9.860 

1.029 


1536 


-3.015 

1. 734 


10.185 

0.966 




-0.563 

5.919 


13.348 

2.872 


1407 


0.992 
1 .293 


9.924 

1.086 




3.053 

5.324 


13. 126 
2.74 5 


1342 


-4.441 
3. 758 


9.464 

3.607 




-0.208 
9. 289 


14.200 
8. 729 


1664 
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TABLE IX 

Keans and Sigmas Beginning Two Hours After Initial Solutions 



J. D. 

61 



68 



69 



Units 



Station 


X mean 
Sigma X 


I mean 
Sigma Y 


Z mean 
Sigma Z 


Point 

Position 


Sample 

Size 


JRK 


6.288 
0.56 1 


-12. 569 
1.571 


3.674 

0.780 


14.597 

1.269 




HM1 


5. 414 
0.492 


-12.937 

1.657 


5.333 

0.689 


15.045 

1.492 


1218 


J HR 


6.29 1 
0. 55 1 


-12. 874 
1. 482 


3.642 

0.770 


14.844 
1. 149 




HH2 


5.315 

0.810 


- 14. 004 
0. 369 


6.578 

1.846 


16.454 

1.056 


1264 


EM 1 


5.407 

0.505 


-13.043 
1. 645 


5.374 

0.688 


15. 150 
1.469 




RM 2 


5. 367 
0. 81 8 


-14. 069 
0. 293 


6.950 
1. 481 


16.650 

0.889 


1150 


BZERO 


5. 51 3 
0. 41 1 


-7. 769 
0. 755 


-2.074 

0.530 


9.780 

0.652 




PILOT 


5. 432 
0.317 


-7. 862 
1. 193 


1.585 

0.764 


9.73 6 
1.072 


1236 


BZERO 


5. 600 
0. 578 


-7. 793 
0.795 


-2.267 

0.739 


9.903 

0.753 




BNOB 


5.476 

0.697 


-1 1.461 
2. 637 


1.997 

3.454 


13.298 

2.802 


1107 


PILOT 


5. 461 
0.337 


-7.910 
1. 293 


1 .552 
0.827 


9 '.795 
1.157 




BNOB 


5.480 
0. 715 


-11.712 
2. 510 


2.414 
3. 114 


13.432 

2.787 


1042 


BZERO 


6.262 
0. 169 


-3. 504 
1. 523 


-3. 053 
2.073 


8.197 

0.501 




PILOT 


5. 717 
0. 940 


-6.001 
4. 678 


3.690 
4. 085 


10.657 

2.850 


1364 



are meters. 
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given X coordinates by approximately 5 to 6 meters at all 
stations on all days of data acquisition having an average 
standard error of .564 meters for all determinations. The Z 
coordinate has smaller discrepancies at all stations on days 
68 and 69 compared to the other coordinates. The average 
standard errors of the Y and Z coordinate discr epancies for 
all determinations are 1.600 and 1.560 meters, respectively 
(Tab. IX). 

The mean point position discrepancies in the KGS-72 
Cartesian coordinate system range between 16.650 meters at 
station RM2 on day 6 1 and 8. 197 meters at station 3ZEE0 on 
day 69. The average of all the mean point position discrep- 
ancy determinations is 12.685 meters. The standard errors 
of the mean point position discrepancies range between .501 
meters at station BZERO on day 69 and 2.802 meters at 
station BNOB on day 68, with an average value of 1.383 
meters for all determinations (Tab. IX) . 

E. G'EODETIC COMPARISONS 

The observed WGS-72 coordinates of each station are 
computed by subtracting the observed discrepancies from the 
given coordinates. The observed WGS-72 Cartesian coordi- 
nates are then transformed to NAD-27 geodetic coordinates 
using a modified version of NGS's D035-41 transformation 
program. 

1 * Discrepancies of Point Position Determinations 

The point position discrepancies of all pairs of 
determinations are computed for the NAD-27 coordinates (Tab. 
X) . Azimuths are measured from north. The average point 
position discrepancy is -.08336 seconds of latitude, -.26088 
seconds of longitude, and *9. 391 meters in ellipsoid height. 
The horizontal position discrepancies of the observed point 
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Geodetic and Horizontal Point Position 
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positions for all determinations range between 6.498 meters 
at station PILOT on day 69 and 8.530 meters at station BZERO 
on day 68. The average distance discrepancy along the 
ellipsoid is 7.596 meters. The azimuths from the given 
positions to the observed positions range between 129° 41’ 

53 ” . 04469 at station BZERO on day 58 and 84° 48' 17". 17181 

at station PILOT on day 69. The observed positions have an 
average azimuth of 109° 25' 56". 5699 from the given 

positions. 

2 • Discrepancies of Point Positions for Entire 

Reco rding Se ssi ons 

There are two determinations of point position 
discrepancies on days 61 and 68 as a result of comparing 
pairs of stations in the computations. An estimate of the 
repeatability of an entire recording session can be calcu- 
lated by first averaging the two point position discrepan- 
cies for each station then computing the standard error of 
the eight determinations. 

The average discrepancies of entire recording 
sessions are -.07980 seconds of latitude and -.25903 seconds 
of longitude with standard errors of .03976 and .01764 
seconds, respectively. The average height discrepancy is 
8.658 meters above the station with a standard error of 
4.748 meters. The average horizontal distance discrepancy 
is 7.538 meters with a standard error of .714 meters. 

3. Relative Pos ition Compa ri sons 

The distances and azimuths between pairs of stations 
are computed for both the given coordinates and the observed 
coordinates by the inverse computation. The differences of 
two are computed by subtracting the given distances and 
azimuths from the observed values (Tab. XI) . Six of the 
seven observed side-lengths are shorter than the distances 
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TABLE XI 



Discrepancies of 


Observed Distances and Azimuths 


J.D. Side 


Dist . (m. ) 


Azimuth 


61 






JRM-RM 1 


1.397 


-44' 44". 2009 


JHR-RM2 


-0.006 


-40 17 * 40". 3464 


RH 1-RM2 


-0.231 


-10 14' 30". 5102 


58 






BZERO-PILOT 


-3. 128 


-1". 24259 


BZERO-BNOB 


-0.428 


-26" .09039 


BNOB-PILOT 


-0.774 


9". 49565 
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BZERO-PILOT 


-4. 556 


-1". 05927 


computed from the given 


coordinates. 


In the short base 


tests, the differences 


in azimuth are 


of the order of 



degrees compared to differences of the order of seconds of 
arc in the medium base line tests. In general, the azimuth 
differences decrease as the side-length differences increase 
in both the short base line and medium base line tests. 

4 . Relative Pos ition Discrepancies 

Relative positioning is used to determine an unknown 
position relative to the position of a known station. The 
direct computation calculates the coordinates of a position 
given the coordinates of a known position and the distance 
and azimuth to the unknown position. The distance and 
azimuth between pairs of observed positions are used in the 
calculation. Stations JMR and RH 1 are taken as fixed in the 
short base line tests and the positions of RM1 and RM2 are 
calculated. Stations BZERO and BNOB are used to locate 
stations BNOB and PILOT in the medium base line tests. The 
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discrepancies of the calculated positions are computed in 
terms of latitude, longitude, elevation, distance, azimuth, 
and relative accuracy between stations (Tab. XII). The 
elevation discrepancies are determined by computing the 
difference of the height discrepancies between pairs of 
point position solutions. 

Relative position determinations differ from the 
given positions by an average of 1.452 meters ir. the short 
base line tests and 2.712 meters in in the medium base line 
tests. The relative accuracies between stations range 
between 1:13 for the side JMR-RM1 on day 61 and 1:17,418 for 
the side BN0B-PI10T on day 68. The average relative accu- 
racy of the medium base line tests is 1:9,712. Station 
PILOT was determined to relative accuracies of 1:8,139 and 
1:5,593 with respect to station ENQ3 on subsequent days. 

The data acquired in the short base line tests 
provide useful information in the absolute position anal- 
yses; however, this scenario is not a practical application 
of relative position determinations. It is not justifiable 
to measure distances of the order of thirty meters with such 
expensive equipment when the distances can be measured accu- 
rately with steel tape or other inexpensive equipment. 
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in Relative Position Determinations 
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V. CONCLUSIONS AND RECOMMENDATIONS 



Based on the results of the medium base line tests, the 
real-time solutions of the TI4 100/3EOSTAR are capable of 
establishing relative positions to Third-order, Class II 
standards. Two receivers, one located at a known position 
and the other separated by 14 to 25 kilometers, can locate 
horizontal control in real-time. This assumes five hours of 
visibility of four satellites and a means for the output of 
a solution based on the averaging of real-time solutions 
recorded at 12-second intervals. 

It is found that relative positioning using GPS is more 
accurate than absolute point positioning. The average 
distance discrepancy of the absolute position determinations 
is 7.538 meters compared to relative position averages of 
1.452 meters in the short base line tests and 2.712 meters 
in the medium base line tests. The absolute point position 
determinations have a one sigma repeatability of .714 meters 
in distance discrepancy. The relatively high repeatability 
of absolute point position determinations resulted in rela- 
tive position accuracies of 1:9,712, on average, for 
distances of 14 to 26 kilometers. 

Third-order, Class II accuracies can be obtained with 
the GEOSTAR Feceiver in real-time. A typical operational 
unit consists of two, 2-man teams; each team supplied with 
the GEOSTAR field set. One team sets an antenna over a 
known position and the other team locates at the point to be 
surveyed. Both receivers should gather data simultaneously 
during the periods of satellite visibility (during Phase II 
testing) . The teams should agree to begin the averaging 
algorithm at a predetermined time tag, preferably a short 
time after the uploads to the satellites. Satellite upload 
times can be requested beforehand from Vandenburg AFB. 
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Additional averaging methods should be explo 
regard to the real-time solutions, averaging ana 
be performed to determine the amount of data acq 
time required to achieve certain accuracies. An 
this nature would be important for special milit 
tions where an accurate position is required in 
five hours of receiver operation. 

In the Fall 1984 the GEOSTAR will undergo ha 
software modification before full-scale producti 
planned that the software will be modified to in 
time solutions of short time-interval observatio 
solutions will be different from existing reai-t 
tions based on a cumulative least squares method 
Relative position determinations using the real- 
tions should improve with the new software becau 
tions of two receivers can be compared simultane 
no influence from preceding solutions. 
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APPENDIX A 

GPS USER-TAPE READ-WRITE PROGRAM 



c 

C THIS PROGRAM IS A SLIGHTLY MODIFIED VERSION OF THE 

C EXAMPLE INPUT PROGRAM CONTAINED IN "A STANDARDIZED 
c EXCHANGE FORMAT FOR NAVSTAR GPS GEODETIC DAI A" BY 

C BY V. DAN SCOTT AND J. GARRY PETERS, APPLIED RESEARCH 

C LABORATORIES, THE UNIVERSITY OF TEXAS AT AUSTIN, MARCH 
C 10, 1983. THIS PROGRAM CAN BE EASILY CHANGED TO READ 

C AND WRITE ANY DATA TYPE ON GPS USER-TAPES SUPPLIED BY 

C NSWC . TO WRITE ALTERNATIVE DATA TYPES, CONSULT THE 
C AEOVE PUBLICATION TO DETERMINE THE DATA TYPES THAT 
C CORRESPOND TO THE DATA OF INTEREST. EXAMPLE: TO PRINT 

C OUT EPHEMERIS INFORMATION (DATA TYPE 4110), CHANGE 

C THE GO TO STATEMENT FOLLOWING THE IF STATEMENT THAT 

C TESTS FOR THIS DATA TYPE BELOW TO TRANSFER CONTROL 

C TO WRITE STATEMENTS WITH THE PROPER FORMAT. THE 

C FORMATS OF ALL DATA TYPES ARE GIVEN IN THE ABOVE 

C PUBLICATION. THE JCL FOR THIS PROGRAM IS FOR AN 
C IBM 370 MODEL 3033. THE TAPES ARE RECORDED WITH 1600 
C BPI, A LOGICAL RECORD LENGTH OF 80, A BLOCKSIZE OF 
C 6400, AND ASCII CODED. 

C 

C*** ** ***** £*$ * *** ** *** **** *** **** ******** * ***** * £ ***** * 

c 

//RAKOW JOB (1643,0130) , ’RAKOWSKY* ,CLASS = F 

//* M A IN ORG = NPGVM1. 1643P 

// EXEC FORTXCLG 

//FORT.SYSIN DD * 

c 

c 
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REAL*8 TTAG,X ,Y,Z,SIGT 



c 

INTEGER IC N , IE EP, ICONI R, NR EP, IN DEX , DTYPE, JUNK , ASTER , 

* BLANK, SENT, PL US, CHARE ,CHARF, CHARI, N RID 
REAL SIGX,SIGY,SIGZ,SKIP 

c 

DATA ASTER/’ *’/, PLUS/’ + ’/, CHARE/’ E’ /, CHAR F/’ F’/, 

* CHAR/’ I’ /, BLANK/’ ’ /, SENT/’ $ ' / 
c 

DIMENSION ICON TR ( 1 9) 

C 

100 CONTINUE 

READ (1,10) ICN,IREP,ICONTR 
10 FOFMAT (A 1, 13, 1 SI4) 

C 

120 CONTINUE 

IF (ICN.EQ. ASTER) GO TO 150 
IF (ICN.EQ. PLUS) GO TO 100 
IF (ICN. EQ. CHARE) GO TO 100 
IF (ICN.EQ. CHARF) GO TO 500 
IF (ICN. EQ. CHARI) GO TO 500 
IF (ICN. EQ. SENT) GO TO 500 
C 

150 CONTINUE 
C 

IF (IBEP. EQ.O) IREP= 1 
NREP=0 
ICN= BLANK 
C 

200 CONTINUE 
C 

INDEX=0 
NREP=NREP+ 1 

IF (IREP.LT.O) GO TO 300 
IF (NREP. GT.IREP) GO TO 100 
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GO TO 300 



C 

300 CONTINUE 
C 

IF (ICH.NE. BLANK) GO TO 120 
I NDE X=INDEX+ 1 
IF (INDEX. GT. 19 ) GO TO 200 
DTYPE=ICONTR (INDEX) 

IF (DTYPE. EQ. 0) GO TO 200 
C 

IF (DTYPE. EQ. 9 1 0) GO TO 204 
IF (DTYPE. EQ. 920) GO TO 201 



IF (DTYPE. EQ. 10 10) 


GO 


TO 


204 


IF (DTYPE. EQ. 10 20) 


GO 


TO 


205 


IF (DTYPE. EQ. 1030) 


GO 


TO 


204 


IF (DTYPE. EQ. 1040) 


GO 


TO 


20 3 


IF (DTYPE. EQ. 10 50) 


GO 


TO 


204 


IF (DTYPE. EQ. 1110) 


GO 


TO 


204 


IF (DTYPE. EQ. 1150) 


GO 


TO 


202 


IF (DTYPE. EQ. 20 10) 


GO 


TO 


20 1 


IF (DTYPE. EQ. 2020) 


GO 


TO 


201 


IF (DTYPE. EQ. 2030) 


GO 


TO 


20 1 


IF (DTYPE. EQ. 21 10) 


GO 


TO 


203 


IF (DTYPE. EQ. 21 20) 


GO 


TO 


20 1 


IF (DTYPE. EQ. 21 30) 


GO 


TO 


201 


IF (DTYPE. EQ. 21 40) 


GO 


TO 


20 1 


IF (DTYPE. EQ. 22 10) 


GO 


TO 


202 


IF (DTYPE. EQ. 2250) 


GO 


TO 


20 1 


IF (DTYPE. EQ. 2270) 


GO 


TO 


201 


IF (DTYPE. EQ. 30 10) 


GO 


TO 


202 


IF (DTYPE. EQ. 3210) 


GO 


TO 


20 1 


IF (DTYPE. EQ. 3220) 


GO 


TO 


201 


IF (DTYPE. EQ. 3250) 


GO 


TO 


201 


IF (DTYPE. EQ. 40 10) 


GO 


TO 


201 


IF (DTYPE. EQ. 4050) 


GO 


TO 


20 1 
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IF (DTYPE. EQ. 4 1 10) 


GO 


TO 


204 


IF (DTYPE. EQ. 4120) 


GO 


TO 


20 2 


IF (DTYPE. EQ. 42 10) 


GO 


TO 


203 


IF (DTYPE. EQ. 4250) 


GO 


TO 


202 


IF (DTYPE. EQ. 4270) 


GO 


TO 


202 


IF (DTYPE. EQ. 50 10) 


GO 


TO 


20 1 


IF (DTYPE. EQ. 60 10) 


GO 


TO 


6010 


IF (DTYPE. EQ. 60 20) 


GO 


TO 


202 


IF (DTYPE . EQ. 6 1 00) 


GO 


TO 


201 


IF (DTYPE. EQ. 6 1 10) 


GO 


TO 


20 3 


IF (DTYPE. EQ. 70 10) 


GO 


TO 


201 


IF (DTYPE. EQ. 71 10) 


GO 


TO 


202 


IF (DTYPE. EQ. 71 20) 


GO 


TO 


201 


IF (DTYPE. EQ. 72 10) 


GO 


TO 


20 1 


IF (DTYPE. EQ. 7220) 


GO 


TO 


201 


IF (DTYPE. EQ. 7230) 


GO 


TO 


201 


IF (DTYPE. EQ. 7240) 


GO 


TO 


202 


IF (DTYPE. EQ. 7250) 


GO 


TO 


20 2 


IF (DTYPE. EQ. 7260) 


GO 


TO 


203 



WRITE (4,31 ) 

31 FORMAT (’ DATATYPE ERROR’) 

C 

201 CONTINUE 

READ (1,11) ICN 
11 FORMAT (A1) 

IF ( (ICN. ME. BLANK) . AND. (ICN. NE. PLUS) ) GO TO 120 
GO TO 300 
C 

202 CONTINUE 

RE AD (1,11) ICN 

IF ( (ICN. NE. BLANK) .AND. (ICN. NE. PLUS) ) GO TO 120 
READ (1,11) NCN 
GO TO 300 
C 

203 CONTINUE 
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READ (1,11) ICN 

IF ( (ICN. NE. BLANK) . AND. (ICN. NE. PLUS) ) GO 10 120 
READ ( 1 , 1 1 ) NCN 
READ (1,11) NCN 
GO TO 300 
C 

204 CONTINUE 

READ (1,11) ICN 

IF ((ICN.NE. BLANK) .AND. (ICN.NE. PLUS) ) GO TO 120 
READ (1,11) NCN 
READ (1,11) NCN 
RE AD ( 1 , 1 1 ) NCN 
GO TO 300 
C 

205 CONTINUE 

READ (1,11) ICN 

IF ( (ICN. NE. BLANK) . AND. (ICN. NE. PLUS) ) GO ID 120 
READ (1,1 1) NCN 
READ (1 ,1 1 ) NCN 
R EAD (1,11) NCN 
READ (1,11) NCN 
GO TO 300 
C 

6010 CONTINUE 



READ ( 1 , 1 2) ICN,NRID,TTAG,X,Y 

IF ( (ICN. NE. BLANK) . AND. (ICN. NE. PLUS) ) GO ID 120 
READ (1 , 13) NCN,Z,SIGT, SIGX,SIGY,SIGZ, SKIP 
WRITE (4, 14) TTAG, X, Y, Z 
WRITE (4 , 15) SIGX,SIGY,SIGZ 
C 

C READ FORMATS 
C 



12 FORMAT (A1 ,12 ,8X,D22. 15,D22. 15, D 17. 10) 

13 FORMAT (A1,2D17.10,4E10. 3) 
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c 

c write foemais 

c 

14 FORMAT (D22. 15, D22. 15, D17. 10 f D 17. 10) 

15 FORMAT (31X,D17. 10,4X,E10.3,4X,E10.3) 

C 

GO TO 300 
500 CONTINUE 
STOP 
END 
c 
c 

//GO. FT0 1F001 DD UNIT=3400-3 , VOL=SER=my tape , 

// DISP= (OLD, PASS) , LABEL= (2, NL, , IN) , 

// DCB= (RECFM = FB, LR ECL=80 , BLKSIZE = 6 400,DEN=3,DPTCD=Q) 
//GO. FT04F001 DD DSN=MSS. SI 643 . FILENAME, UNIT=333 0V , 

// DISP= (NEW, CATLG, DELETE) , MSVGP=PUB4B, 

// DCB= (RECFM = FB, LR ECL=80 , BLKSIZE = 64 00 ) 

/* 

// 
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APPENDIX 3 

PROGRAM TO COMPOTE REAL-TIME SOLOTION DISCREPANCIES 



£**********#********************************************** 

c 

C THIS PROGRAM COMPUTES THE DIFFERENCES OF THE X, Y, 

C AND Z COORDINATES, THE TOTAL POINT POSITION DISCREP- 
C ANCIES, AND THE SIDE-LENGTH DISCREPANCIES FOR TWO 
C STATIONS AT COMMON GPS TIMES. THIS PROGRAM, IN EFFECT, 
C WRITES THE ABOVE DATA BEGINNING WITH THE FIRST REAL- 
C TIME SOLUTIONS OF EACH TAPE, THEN SPOOLS THE TAPE OF 
C THE EARLIER RECORDING SESSION TO CORRESPOND TO THE 

C GPS TIME OF THE TAPE THAT STARTED RECORDING AT A LATER 

C TIME. A CONVERSION OF THE GIVEN GEODETIC COORDINATES 
C TO V7GS-72 CARTESIAN COORDINATES MUST FIRST BE DONE. 

C THESE VALUES ALONG WITH COMPUTED SIDE-LENGTH DISTANCE 
C ARE ENTERED BELOW FOLLOWING THE TYPE DEFINITIONS. THE 
C OUTPUT IS IN THE FORM: ’GIVEN VALUE - REAL-TIME 

C SOLUTION VALUE’ . 

C 

C*** ’S'********'***** ***** ** ** ******************************* 
C 

//RAKOW JOB (1643, 01 30} , ’ RAKOWSKY’ ,CLASS=B 
//* MAIN 0RG=NPGVM1. 1643P 
// EXEC FORTXCLG 
//FORT. SYSIN DD * 
c 

REAL*8 TTAG1 , TTAG2, XI ,X2,Y1,Y2,Z1,Z2, DSL X , DEL Y , D ELZ , 

* DIST,D VALUE, SIGX1 ,SIGX2,SIGY1 ,SIGY2, SI3Z 1, SIGZ2, 

* DELDST, X 1GIVN,X2GIVN, Y1GIVN, Y23I VN , Z 1 GI V N , Z2GI V N, 

* XI DIF, X2DIF, Y1DIF, Y2 DI F, Z 1 DI F, Z 2 DIF, 0 FSET 1 , OFSET2 
C 

D VALUE=25 493. 1 88 
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X lGIVN=-740577.548 
YlGIVN=- 545677 3. 662 
Z1GIVN =3207660 .08 1 
X2GIVN=-740336 .50 1 
Y2GIVN=-54697 17.020 
Z2GIVN=3185698.432 
C 

REAL (1,10, END=500) TTAG 1 , X 1 , Y 1 , Z 1 
R EAD (1,11) SIGX1,SIGY1 , SIGZ1 
READ (2,10, END= 500) TTAG2, X2, Y2, Z2 
R EAD (2,11) SIGX2,SIGY2 ,SIGZ2 
C 

IF (TTAG1.LT. TTAG2) GO TO 3 
C 

2 CONTINUE 

R E AD (2 , 10) TTAG2,X2,Y2,Z2 
READ (2,11) SIGX2,SISY2,SIGZ2 
IF (TTAG2.GE. TTAG1) GO TO 30 
GO TO 2 
C 

3 CONTINUE 

READ (1,10) ITA G 1 , X 1 , Y 1 , Z 1 
R EAD (1,11) SIGX1,SIGY1 , SIGZ1 
IF (TTAG1 .GE.TTAG2) GO TO 30 
GO TO 3 
C 

30 CONTINUE 

IF (TTAG1.GT.TTAG2) GO TO 2 

IF (TTAG2.GT. TTAG1) GO TO 3 

DE1X=X1-X2 

DELY=Y 1 - Y2 

DE1Z=Z1-Z2 

DIS7=DS2RT (DELX**2+DELY**2+DELZ**2) 
DELD5T=DIST-DVALUE 
C 
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c 



c 



X1DIF=X1-X1GIVN 
Y 1DIF=Y1-Y1GIVN 
Z1DT T ’=Z1-Z1GIV N 
X2DIF=X2-X2GIVN 
Y2DIF=Y2-Y2GIVN 
Z 2DI F=Z2- Z 2GI VN 

0FSET1 =DSQRT (X1DIF**2 + Y1DIF**2+Z1DIF**2) 
0FSET2=D SQ RT (X2DIF**2 + Y2DIF**2 + Z2DIF**2) 



WRITE (4, 42) 
WRITE (4,43) 
WRITE (4, 4 1) 
WRITE (4,43) 
WRITE (4,4 1) 
WRITE (4,41) 



TTAG1 , TTAG2 

XI DIF,Y1DIF, ZlDIF 

0FSET1 

X2DIF, Y2DIF, Z2DIF 

0FSET2 

DELDST 



READ (1 , 1 0, END=500) TTAG1,X1,Y1,Z1 
R EAD (1,1 1 ) SIGX1,SIGY1 ,SIGZ1 
READ (2,1 0, END=500) TTAG2,X2,Y2,Z2 
R EAD (2,11) SIGX2,SIGY2,SIGZ2 
GO TO 30 



C 

C READ FORMATS 
C 



10 FORMAT (D22. 15,D22. 15,D17. 10, D 17. 10) 

11 FORMAT (31X,D17 . 10,4X,E10.3,4X,E10. 3) 
C 

C WRITE FORMATS 
C 

41 FORMAT (D22. 15) 

42 FORMAT (2D22. 15) 

43 FORMAT (3D22. 15) 

C 

500 CONTINUE 
STOP 
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EIT D 



C 

//30. FT01F001 DD DSN=MSS. S1643. MYTAPE, 

// DISP= (OLD, KEEP) ,MSVGP=PUB4S 

//SO. FT02F001 DD DSN=MSS. SI 643 .FILENAME, 

// DISP= (OLD, KEEP) ,MSVGP=PU34B 

/*30. FT04F001 DD D SN =MSS - SI 643 . CHEZ K , UNI T=3 3 50 , 

/* DISP= (NEK, KEEP) ,VOL=SER=MVS004,SPACE= (CYL, (4,4)) , 

/* DC E= (RECFM =FB, LRECL=80 , BLKSI ZE = 640 0) 

/* 

//SO. FT04F001 DD DSN=MSS. SI 643. filename, UNIT=3330V , 

// DISP= (NEW, CATLG, DELETE) , MS VGP = PU34 B, 

// DCB= (EECFM = FB,LEECL = 80,BLKSIZE = 6400) 

/* 

// 
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APPENDIX C 

LIST OF ABBREVIATIONS 



AODC 


Age of Data (Clock) 


AOOE 


Age of Data (Ephemeris) 


ARL:UT 


Applied Research Laboratories: Univ. of Texas 


ASCII 


American Standard Code for Info. Interchange 


CD(J 


Control and Display Unit 


DMA 


Defense Mapping Agency 


DMAHTC 


Defense Mapping Agency Hydro. /Topo. Center 


GEOSTAR 


Geodetic Satellite Tracking and Ranging 


GMT 


Greenwich Mean Time 


GPS 


Global Positioning System 


Mbps 


Million bits per second 


MCS 


GPS Master Control Station 


NAD-27 


North American Datum of 1927 


NAVSTAR 


Navigation Satellite Tracking and Ranging 


NSWC 


Naval Surface Weapons Center 


NSWC/DL 


Naval Surface Weapons Centar/Dahlgren Labs. 


PPS 


Precise Positioning Service 


PR N 


Pseudo-Random Noise 


SPS 


Standard Positioning Service 


WGS-72 


World Geodetic System of 1972 
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